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ABSTRACT 

We present rotational and astrometric parameters of three millisecond pulsars located near the center of the 
globular cluster NGC 6266 (M62) resulting from timing observations with the Parkes radio telescope. Their 
accelerations toward the cluster center yield values of the cluster central density and mass-to-light ratio consistent 
with those derived from optical data. The three pulsars are in binary systems. One (spin period P =5.24 ms) is in 
a 3.5-day orbit around a companion of minimum mass 0.2 Mq. The other two millisecond pulsars (P =3.59 ms 
and 3.81 ms) have shorter orbital periods (3.4 h and 5.0 h) and lighter companions (minimum mass 0.12 Mq and 
0.07 Mq respectively). The pulsar in the closest system is the fifth member of an emerging class of millisecond 
pulsars displaying irregular radio eclipses and having a relatively massive companion. This system is a good 
candidate for optical identification of the companion star The lack of known isolated pulsars in NGC 6266 is also 
discussed. 

Subject headings: Globular clusters: individual (NGC 6266) — pulsars: individual (PSR J1701-3006A, 
PSR J1701-3006B, PSR J1701-3006C) 



1. INTRODUCTION 

Recycled pulsars are old neutron stars revived through trans- 
fer of matter and angular momentum from a mass-donor com- 
panion in a binary system (e.g. Alpar et al. 1982; Smarr & 
Blandford 1976; Bhattacharya & van den Heuvel 1991). They 
are point-like objects and can be considered as test masses for 
probing gravitational effects. Most of them are also extremely 
stable clocks, allowing for accurate measurements of their rota- 
tional parameters, position and apparent motion in the sky. Be- 
cause of these characteristics, recycled pulsars found in globu- 
lar clusters (GCs) have proven to be valuable tools for studying 
the GC potential well (e.g. Phinney 1992; Camilo et al. 2000; 
D'Amico et al. 2002), the dynamical interactions in GC cores 
(e.g. Phinney & Sigurdsson 1991; Colpi, Possenti & Gualan- 
dris 2002) and neutron star retention in GCs (e.g. Rappaport et 
al. 2001). In the case of 47 Tucanae they allowed also the first 
detection of gas in a GC (Freire et al. 2001). 

Globular clusters are a fertile environment for the formation 
of recycled pulsars: besides evolution from a primordial sys- 
tem, exchange interactions in the ultra-dense core of the cluster 
favor the formation of various types of binary systems suitable 
for spinning up the neutron stars they host (Davies & Hansen 
1998). Because of this, about 60% of all known millisecond 
pulsars (MSPs) are in GCs. Unfortunately, pulsars in GCs are 
elusive sources since they are often distant and in close binary 
systems. Their large distances make their flux density typically 
very small and their signals strongly distorted by propagation 
through the dispersive interstellar medium. In addition, they 



frequently are members of close binary systems, causing large 
changes in the apparent spin period and sometimes periodic 
eclipsing of the radio signal. 

The Parkes Globular Cluster (PKSGC) survey is a search for 
pulsars in the system of southern GCs using the Parkes 64-m 
radio telescope which commenced in 2000. It exploits the high 
sensitivity of the central beam of the Parkes multibeam receiver 
(Staveley-Smith et al. 1996), the efficiency of a modern data 
acquisition system (e.g., Manchester et al. 2001, D'Amico et 
al. 2001a) and the high resolution of a new filterbank designed 
and assembled at Jodrell Bank and Bologna, with the aim of 
improving the capability for probing distant clusters. Time se- 
ries data are analyzed with a modern algorithm for the incoher- 
ent search of periodicities over a range of dispersion measures 
(DMs) and accelerations resulting from orbital motion. 

This project has akeady resulted in the discovery of 12 mil- 
lisecond pulsars in six globular clusters which had no previ- 
ously associated pulsar (D'Amico et al. 2001a; D'Amico et al. 
2001b; Possenti et al. 2001; D'Amico et al. 2002). These de- 
tections reversed the declining trend in discoveries of additional 
clusters hosting these objects; in the seven years from 1987 
(when the first pulsar in a GC, B 182 1-24 in M28, was discov- 
ered at Jodrell Bank by Lyne et al. 1987) to 1994 (B1820-30A 
and B in NGC 6624: Biggs et al. 1994) 13 globular clusters 
were shown to contain at least one pulsar, whereas no new clus- 
ter joined the list in the following six years. More recently, pul- 
sars have been detected in a further three GCs (Ransom 2003a; 
Ransom 2003b; Jacoby 2003), bringing the current total to 73 
pulsars in 22 clusters.^ 



The association of the long-period pulsar B1718— 19 with the cluster NGC 6342, questioned by some, is included in this list. 
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This paper discusses results from the PKSGC survey of 
NGC 6266 (M62). The discovery of the first pulsar in this 
cluster, PSR J1701-3006A, was presented by D'Amico et al. 
(2001a). A preUminary announcement of the discovery of 
two more MSPs was also made by D'Amico et al. (2001b), 
while three further millisecond pulsars were later detected at 
the Green Bank Telescope (Jacoby et al. 2002). Here we re- 
port details of the discovery of the second and third pulsars, 
PSRs J1701-3006B and J1701-3006C, both members of short- 
period binary systems, and discuss timing results obtained over 
a 3-yr interval for all three MSPs discovered at Parkes. Based 
on these results we investigate the properties of the pulsars and 
of the host cluster. We particularly discuss the MSP in the 
tightest of the three systems, which belongs to the rare class 
of ecUpsing radio pulsars. 

2. DATA COLLECTION AND PROCESSING 

The PKSGC survey uses the dual polarization central beam 
of the 20-cm multibeam receiver of the Parkes radio telescope. 
The two channels have a system temperature of ~ 22 K and a 
central frequency of 1390 MHz. A high-resolution filterbank 
system consisting of 512 x 0.5 MHz adjacent channels per po- 
larization is used to minimize dispersion smearing, preserving 
significant sensitivity to a 3 ms pulsar with dispersion measure 
up to 300 cm~^pc. In this case, the limiting sensitivity is ~ 0.15 
mJy for a signal-to-noise ratio (s/n) of 8 and a standard 2-h 
observation (assuming a typical duty cycle of ~ 20% and negli- 
gible scattering). After adding the outputs in polarization pairs, 
the resulting 512 data streams are each high-pass-filtered, in- 
tegrated and 1-bit digitized every 125 /zs. Each observation 
typically produces 2-4 Gbytes of data; a cluster of 10 Alpha- 
500MHz CPUs at the Astronomical Observatory of Bologna 
has been used for offline processing. 

The processing pipeline first splits each data stream into non- 
overlapping segments of 1050, 2100, 4200 or 8400 s, which are 
processed separately. When no pulsar is known in a GC (and 
the DM is therefore unknown) the data are dedispersed over a 
wide range of ~ 500- 1000 trial DMs, spanning the interval 
(1.0 ±0.4) DMexp, where DMexp is the DM expected for the 
cluster according to a model for the Galactic distribution of the 
ionized gas (Taylor & Cordes 1993). Each dedispersed series 
is then transformed using a Fast Fourier Transform, and the re- 
sulting spectra are searched for significant peaks. The process 
is repeated for spectra obtained from sunaming 2, 4, 8 and 16 
harmonics. This produces a large number of candidate periods 
above a threshold. The time-domain data are then folded in 
sub-integrations at each of these periods in turn and searched 
for both a linear and a parabolic shift in pulse phase. A hn- 
ear shift corresponds to a correction in the candidate period, 
whereas a parabolic correction is a signature of acceleration of 
the pulsar due to its orbital motion. Parameters for final pulse 
profiles with significant s/n are displayed for visual inspection. 
This processing scheme led to the discovery of the first pulsar 
in NGC 6266 (D'Amico et al. 2001a). 

For an MSP in a very close orbit and with relatively high 
minimum companion mass, the acceleration may undergo sig- 
nificant variations during an observation. As a consequence, 
weak sources can be missed at the confirmation stage if a con- 
stant acceleration term is applied to the data. Hence code has 
been developed at Bologna for searching both the acceleration 
and the derivative of the acceleration in the sub-integration ar- 

^See http://www.atnf.csiro.au/research/pulsar/timing/tempo. 



rays of interesting candidates. Spanning a cubic space, the code 
searches also for the period in a small interval of values around 
the nominal candidate period. Using this code we were able to 
confirm two more binary millisecond pulsars in NGC 6266. 

Once a pulsar is detected and confirmed in a cluster, the data 
are reprocessed with dedispersion at the DM value of the newly 
discovered pulsar The resulting time series is then subject 
to a fully coherent search for Doppler-distorted signals over 
a large range of acceleration values. Applying this extremely 
CPU-intensive procedure to NGC 6266, exploring accelerations 
in the interval |a| < 35 m s"^ for 35-min long segments (and 
\a\ < 17.5 m s'^ for 70 -min long segments), resulted in no fur- 
ther discoveries. 
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Fig. 1 . — Post-fit timing residuals as a function of the Modified Julian Day of 
observation (left) and integrated pulse profiles at a central frequency of 1390 
MHz (right) for the three rmllisecond pulsars in NGC 6266 which are dis- 
cussed in this paper. The short horizontal line on the left side of each pulse 
profile represents the time resolution of the integrated profile including DM 
smearing. 

Regular pulsar timing observations at the Parkes 64-m ra- 
dio telescope began shortly after the discovery of these pulsars, 
using the same observing system as the search observations. 
Timing observations, typically of 30 to 60 minutes duration, are 
dedispersed and synchronously folded at the predicted topocen- 
tric pulsar spin period in an off-Une process, forming pulse 
profiles every few minutes of integration. Topocentric pulse 
times of arrival (TOAs) are determined by convolving a stan- 
dard high s/n pulse template with the observed pulse profiles 
and then analyzed using the program tempo^. tempo con- 
verts the topocentric TOAs to solar-system barycentric TOAs 
at infinite frequency (using the DE200 solar-system ephemeris, 
Standish 1982) and then performs a multi-parameter fit to de- 
termine the pulsar parameters. 

Table 1 lists the timing parameters obtained for the three pul- 
sars, including precise positions. Values of the dispersion mea- 
sure (DM) were obtained for each pulsar by splitting the to- 
tal bandwidth into four adjacent 64-MHz wide sub-bands and 
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computing the differential delays. The available data do not 
yet allow a constraining fit for the orbital eccentricity e for 
any of the three pulsars (see the footnotes to Table 1 for de- 
tails of the fitting procedure.) The mean flux densities at 1400 
MHz (5i40()) in Table 1 are average values, derived from the 
system sensitivity estimate and the observed s/n. In the case 
of PSR J1701-3006B, the quoted flux density refers only to 
epochs away from the ecUpse (see below). As expected from 
the relatively high DMs, interstellar scintillation does not sig- 
nificantly affect the detectability of these sources; observed 
variations are within 30% of the nominal flux density reported 
in Table 1. 
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Fig. 2. — Post-fit timing residuals as a function of orbital phase for the three 
milUsecond pulsars in NGC 6266 discussed in this paper. AH the orbits have 
been uniformly sampled, with the exception of PSR J1701-3006B for which 
we have excluded from the fit the TOAs in the region of the eclipse. 



The inferred radio luminosities of the three millisecond pul- 
sars 10-20 mJy kpc^ at 1400 MHz, corresponding to a lu- 
minosity at 400 MHz L400 ;^ 100 mJy kpc'^ for a typical spectral 
index -1 .7, see Table 1 ) places all these sources in the bright tail 
of the luminosity function of millisecond pulsars in the Galac- 
tic disk (Lyne et al. 1998) and in 47 Tucanae (Camilo et. al 
2000). If we assume a luminosity distribution dJV oc L"MlogL 
(Lorimer 2001), NGC 6266 would contain a few hundred pul- 
sars with L400 ~ 1 mJy kpc^, the approximate limiting lumi- 
nosity observed for Galactic disk pulsars. Unfortunately, the 
cluster distance and the lack of any strong signal enhancement 
due to scintillation will make difflcult detecting the fainter pul- 
sar population, probably preventing a direct investigation of the 
shape of the pulsar luminosity function in this cluster. 

The pulsar PSR n701-3006A has flie largest flux density 
and the longest orbital period among the three and was first de- 
tected in an observation during 1999 December (D'Amico et 
al. 2001a). PSR J1701-3006B and J1701-3006C are weaker 
pulsars in closer binary systems and were confirmed in 2000 
November. Once the orbits were determined, signals from these 
two pulsars were recovered in all the observations performed 



prior to their confirmation. Therefore the timing solutions re- 
ported in Table 1 (and whose residuals are displayed in Fig. 1) 
cover the same time-span for all three MSPs, from 2000 June 
to 2003 May. Inspection of Figure 2 shows that all the orbits 
have been uniformly sampled (excepting PSR J1701-3006B for 
which we have excluded TOAs in the region of the eclipse; see 
below) and that there are no systematic trends in the residuals 
as a function of binary phase. 

3. CONSTRAINTS ON PULSARS AND CLUSTER PARAMETERS 

NGC 6266 is listed in flie Webbink (1985) catalog as a 
moderately reddened, E(B-V) = 0.48, medium-low metallic- 
ity globular cluster, with [Fe/H] = -1.38 ±0.15, located at 
~ 6.9 ± 1 .0 kpc from the Sun (Brocato et al. 1996) and proba- 
bly having a collapsed core (Harris 1996). 

The three millisecond pulsars discussed here are all located 
close to the center of mass of the cluster, at least in projec- 
tion, with projected distances < 1.8 8e, where = 10"8 is the 
core radius of NGC 6266 (Harris 1996). This is consistent with 
the hypothesis that the cluster has reached thermal equilibrium, 
in which energy equipartition gives less velocity to the most 
massive species, constraining them to reside deep in the cluster 
potential well. 

The spin period derivatives P are all negative, implying that 
the line-of-sight acceleration a/ imparted to the pulsars is di- 
rected toward the observer and that it overcomes the (positive) 
Pi due to intrinsic spin-down (see e.g. Phinney 1993). The 
probability that a nearby passing star in the crowded cluster 
core is significantly accelerating at least one of the three MSPs 
is < 1% (Phinney 1993). Moreover, given the position and the 
kinematics of the globular cluster NGC 6266, the centrifugal 
acceleration of the system (Shklovskii 1970) and the vertical ac- 
celeration in the Galactic potential (Kuijken & Gilmore 1989) 
produce only negligible effects on the measured a/ = \cP/P\. 
The differential Galactic rotation (Damour & Taylor 1991) can 
contribute at most a positive ^10%, ^2% and ^25% to the ob- 
served P/P of the pulsars A, B and C respectively; hence, we 
conclude that the sign of the line-of-sight accelerations is dom- 
inated by the radiaUy synmfietric mean gravitational field of the 
globular cluster and that the three MSPs are located behind the 
plane of the sky through the cluster center 

The maximum possible o; due to the mean gravitational field 
in NGC 6266 is given by the following relation (accurate at the 
10% level for Gx ^ 20)^, Phinney 1992) 



(3/2) 



D (62 + 92)1/2 



(1) 



where ct; = 14. 3 ±0.4 km s"' is the line-of-sight velocity disper- 
sion (Dubath et al. 1997) and Z) = 6.9 ± 1 .0 kpc is the distance 
(Brocato et al. 1996). @c and 0_l are the angular core radius 
and the angular displacement with respect to the globular clus- 
ter center, located at R.A. (J2000): 17i>01"12?8, Dec. (J2000): 
-30°06'49" (Harris 1996, catalog revision 2003). In particular, 
for a pulsar with negative P the following inequality must hold: 
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where c is the speed of light. 

The observed lower limit on the magnitude of the line-of- 
sight acceleration of PSR J1701-3006B, a, = 2.9 x lO-*" cm s'^, 
is flie fliird largest after those of PSR B2127-I-1 lA andD in M15 
(Anderson et al. 1990) and is almost identical to those of the 
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two MSPs with negative P recently discovered in the central 
regions of NGC 6752 (D'Amico et al. 2002). For NGC 6752, 
the high values of P imply a central mass-to-Ught ratio larger 
than that from optical estimates (D'Amico et al. 2002). For 
NGC 6266 on the other hand, the upper panel in Figure 3 shows 
that the parameters derived from optical observations can en- 
tirely account for the large P/P of PSR J1701-3006B (the ver- 
tical size of the dots in Figure 3 represents the contribution to 
ai/c due to the differential Galactic rotation). In particular, ap- 
plying equation (1) of D'Amico et al. (2002) we derive a lower 
limit on the central mass-to-light ratio (expressed in solar units) 
for NGC 6266, M./L=\ .6, which is compatible with the opti- 
cal value reported in the literature, 2.0 (Pryor & Meylan 1993). 
Similarly, using the observed P/P of J1701-3006A (corrected 
for the Galactic contribution) and equation (7) of Camilo et al. 
(2000) the inferred lower limit po = 2.1 x 10^ MqPC"^ of the 
central mass density of NGC 6266 is within the limits obtained 
from optical data (Pryor & Meylan 1993). These results suggest 
that, even though all three clusters display a compact core and 
very high line-of-sight accelerations for the embedded pulsars, 
the dynamics in the inner region of NGC 6266 probably more 
resemble those of Ml 5, for which 2 < M./ C <3 was inferred 
by Phiimey (1993). 
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Fig. 3. — Upper panel: maximum line-of-sight acceleration |a/,max/c| = 
\P/P\ versus displacement Q± with respect to the center of NGC 6266. The 
solid and the two dashed lines represent the predictions based on equation (1) 
using the nominal values of the distance and of the line-of-sight dispersion 
velocity and their Icr uncertainties obtained from the available optical obser- 
vations (see text). The dot-dashed vertical line marks the assumed angular core 
radius = 10"8 (Harris 1996). The points represent lower limits to the line- 
of-sight accelerations based on the timing solutions for the three millisecond 
pulsars. The vertical size of the points corresponds to the contribution to \P/P\ 
due to the Galactic potential. Lower panel: constraints on the age of PSR 
J1701-3006B obtained from equation (1) and (2). The thin solid Unes and the 
dotted lines represent the values of the parameters reported in literature and 
their \a uncertainties. An intrinsic characteristic age of PSR J1701-3006B 
larger than about 1 Gyr is compatible with the available observations. 

The satisfactory match between the dynamical parameters of 
NGC 6266 constrained from pulsar timing observations and de- 



rived from optical data allows use of the latter for deriving re- 
liable constraints on the age and surface magnetic field of the 
milUsecond pulsars. For instance, the lower panel in Figure 3 
shows that the intrinsic characteristic age of PSR J1701-3006B 
should be greater than ^1.3 Gyr to be consistent with the clus- 
ter's distance and velocity dispersion (including their Icr uncer- 
tainties). This in turn implies an upper limit on the surface mag- 
netic field fij = 3.2 X lO^'^iPPyl^ = 4.0 x 10^ G. Less stringent 
limits can be similarly derived for PSR J1701-3006A (r, > 0. 15 
Gyr and B, < 17 x 10* Gauss) and PSR J1701-3006C (r,- > 25 
Myr and ^ 31 x 10** Gauss). These values are typical for 
MSPs, both in GCs and in the Galactic disk. 

4. RANGE OF DISPERSION MEASURES 

The MSPs in NGC 6266 show the second largest range in 
DM (a maximum deviation ADM = 0.9 cm'-'pc with respect 
to the average DM^ve = 1 14.34 cm'^pc) after PSR B1744-24A 
and PSR J1748-2446C in Terzan 5 (Lyne at al. 2000). This 
large range is probably due to a significant gradient in the 
Galactic electron column density across different lines-of-sight 
toward the cluster, an interpretation supported by the strong 
variations in reddening observed across this cluster: (5E = 
AE(B-V)/E(B-V) = 0.19/0.48 for an angular displacement 
of A^E ~ 7' (as derived from Fig. 3 of Minniti, Coyne & Claria 
1992). Alternately, the variations may have a local origin as in 
the case of 47 Tucanae, where pulsar timing observations show 
that they are due to a plasma permeating the cluster (Freire et 
al. 2001). The same explanation has been proposed in the case 
of M15 (Freire et al. 2001) and NGC 6752 (D'Amico et al. 
2002). For NGC 6266, the electron number density of a uni- 
form fully ionized gas would be surprisingly high, = 1 .6 ±0.4 
cm"^, at least an order of magnitude larger than that estimated 
for the other clusters. The determination of positions, acceler- 
ations and precise DMs of the additional MSPs discovered in 
this cluster (Jacoby et al. 2002) will help in determining the 
origin of the scatter in DM. 

5. THE ABSENCE OF ISOLATED PULSARS IN NGC 6266 

In contrast to other globular clusters in which at least five 
pulsars have been discovered (in order of decreasing number of 
pulsars: 47 Tuc, M15 and NGC 6752), all the MSPs known in 
NGC 6266 are in binary systems (including the three detected 
by Jacoby et al. 2002). The absence of known isolated pulsars 
in NGC 6266 cannot simply be ascribed to a selection effect 
since, for a given spin period and flux density, an isolated MSP 
is easier to detect than a binary MSP. Unfortunately, the ob- 
servational biases affecting the fraction .Fis of isolated pulsars 
discovered in a given cluster (with respect to the total observed 
MSP population) are difficult to quantify precisely. Consid- 
ering all the other clusters, J^i^ ^2/5. If this ratio applies to 
NGC 6266, the probability of having the first six detected pul- 
sars be all binary is ^ 5%. 

If this absence of isolated pulsars in NGC 6266 is not a statis- 
tical fluctuation, it must relate to the mechanisms of formation 
of these objects and their interplay with the dynamical state of 
the cluster. The few isolated milUsecond pulsars observed in 
the Galactic disk (where .Fis ~ 1 / 3) are thought to be endpoints 
of a rare process of ablation and eventually evaporation of the 
companion star by the energetic flux of particles and electro- 
magnetic waves emitted by the pulsar (e.g. Ruderman, Shaham 
& Tavani 1989). Besides this formation channel, the isolated 
MSPs seen in globular clusters can also result from close stel- 
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lar encounters disrupting a binary system which had previously 
been through the recycling process (e.g. Sigurdsson & Phinney 
1993). 

This suggests that NGC 6266 is now in a dynamical state 
where the rate T^^form of formation (and of hardening) of binary 
systems containing a neutron star (and suitable for producing 
new MSPs) is much larger than the rate of disruption T^disr of 
such systems. This idea is supported by comparison of the rel- 
evant rates with other clusters. Table 2 summarizes the values 
of 7?.form and 7?.disr for the four clusters containing at least five 
known pulsars. 7?.form scales as the rate of close encounters in 
the cluster, in turn proportional to Po^r^, where po is the central 
luminosity density and the core radius of the cluster (Ver- 
bunt 2003). Inspection of Table 2 shows that the expected fre- 
quency of close encounters in the core of NGC 6266 ^ is 40% 
larger than that of M15 and seven times that of NGC 6752. 
Although the numbers of known pulsars in these clusters are 
similar, the comparison of their radio luminosities (see discus- 
sion in § 2) indicates that NGC 6266 hosts many more pulsars 
than NGC 6752, in accordance with the trend suggested by the 
values of T^om- 

On the other hand, the probability that a binary, once formed, 
will experience a further encounter, which may change or split 
it (sometimes creating an isolated millisecond pulsar), scales as 
"^disr oc Po^r"' (Verbunt 2003). Hence large values of the ratio 

'^form / "^disr 

OC porl should indicate that more neutron stars are 
in binary systems than are isolated. As Table 2 shows, this pre- 
diction roughly conforms with the numbers for the three globu- 
lar clusters having collapsed cores: NGC 6266 (six binary pul- 
sars) has a ratio T^fom/'T^disr a few times larger than M15 (one 
binary and seven single pulsars) and NGC 6752 (one binary and 
four single pulsars). 

Despite these encouraging agreements, the scaling relations 
may miss many factors which could strongly differentiate the 
pulsar population in the clusters, e.g., the mass-to-light ratio, 
the mass function in the core, the neutron-star retention frac- 
tion, the period distribution of the binary systems and the effects 
of the collapse of the core. The last point could be especially 
relevant, as the only non-coUapsed cluster containing more than 
five known pulsars, 47 Tucanae, fits with the predictions based 
on 7?.form: but does not satisfy those related to T^disr (see Ta- 
ble 2) — its binary disruption rate should be less than half that 
of NGC 6266, but it hosts several isolated millisecond pulsars, 
with a value of J^is similar to that seen in the Galactic field, 
where dynamical encounters are unimportant in the formation 
of isolated millisecond pulsars. 

Detailed numerical simulations are required to investigate if 
trapping of almost all the neutron stars in close binary systems 
can really occur, for instance, during the phase immediately 
preceding the core collapse or its reversal. 

6. THE ECLIPSES IN PSR J1701-3006B 

PSR J1701-3006B displays partial or total eclipses of the ra- 
dio signal at 1.4 GHz near superior conjunction, i.e. at orbital 
phase 0.25 (see Fig. 4), clearly due to gas streaming off the 
companion. A typical event starts at orbital phases in the range 
0.15-0.20 and ends at orbital phase ~ 0.35, hence sometimes 
displaying a slight asymmetry with respect to the expected 
nominal center of the eclipse at phase 0.25. At both eclipse 
ingress and egress, the pulses usually exhibit excess propaga- 

'The absence of pulsars in long period orbital systems (easily destroyed in 
significant rate in the central region of NGC 6266 since a time at least comparable 



tion delays (see Figs. 4 and 5). The eclipse region covers up 
to 20% of the entire orbit but, as illustrated in Figure 4, unpre- 
dictable irregularities affect both the duration and the appear- 
ance of the ecUpses. Sometimes the pulsation remains barely 
visible (see e.g. the case of Fig. 4a), while on other occasions 
the pulse is totally eclipsed for a large portion of the event (e.g. 
the cases of Fig. 4e). In a favorable case (Fig. 4b), it has been 
possible to measure a shght reduction of the s/n of the pulse (al- 
though at the IfT level only: see caption of Fig. 5) as the pulsar 
signal crosses the region of interaction with the matter released 
by the companion. 
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Fig. 4. — Observed signal intensity as a function of orbital phase and pul- 
sar phase for five observations of PSR J1701-3006B centered at 1390 MHz 
with a bandwidth of 256 MHz. Eclipses are expected to occur around superior 
conjunction (phase 0.25). The data are processed in contiguous integrations 
of 120 s duration, (a) ~ 68 min observation starting on 2002 November 27 at 
05:41 UT; (b) ~ 210 min observation starting on 2003 January 26 at 00:01 UT; 
(c) ~ 69 min observation starting on 2000 July 21 at 07:54 UT; (d) ~ 68 min 
observation starting on 2002 July 10 at 07:12 UT; (e) ~ 131 min observation 
starting on 2002 April 29 at 13:52 UT. 

Pulse broadening and reduction of the s/n prevent investiga- 
tion of the frequency-dependent behavior of the delays in our 
256-MHz bandwidth. However, assuming that they are com- 
pletely due to dispersion in an ionized gas (as shown for other 
eclipsing pulsars, e.g. Fruchter et al. 1990; D'Amico et al. 
2001c), the corresponding electron column density variations 
ANe may represent a first viable explanation of the eclipse phe- 
nomenology. With ANe 1.5 X 10^^ At-3 cm~^, where At-3 is 
the delay at 1 .4 GHz in ms, whenever At-T, < 2 (which could 
be the case for the entire events in Figs. 4a and 4b), the implied 
pulse broadening over the receiver bandwidth AP-j, = 0.36 A f_3 
ms is at most 80% of the intrinsic pulse width (~ 0.50 P at 
10% of the peak). Hence the pulse may be only largely broad- 
ened (with an implied reduction of s/n), but not disappear com- 
pletely. On other occasions, the delays may increase much 
more rapidly, possibly growing well beyond At-3 = 2. In this 

dynamical interactions) suggests that close encounters have been occurring at a 
to the cluster relaxation time, ~ 1.5 Gyr at the half-mass radius (Harris 1996). 
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case, the DM variations alone could completely smear the sig- 
nal, causing a total disappearance of the pulsations. 

Alternatively, free-free absorption of the radio-waves in an 
ionized envelope of matter released from the companion and 
expanding adiabatically can explain both the weakening and 
the total disappearance of the radio signal. The optical depth 
for this process can be written (see Spitzer 1978 and Rasio, 
Shapiro & Teukolsky 1989) as 



Tff = 0.74 



1.32 Rp 



0.8 Re 
Re 



(3) 



where the orbital separation a and the radius of the eclipse, 
Re, defined to be the chord at radius a subtended by the angle 
between the orbital phase of eclipse ingress and orbital phase 
0.25, are scaled for PSR J1701-3006B (assuming an orbital 
incUnation of 60 degrees, see later), T is the temperature of 
the fuUy ionized gas and A?_3 is the observed delay in mil- 
liseconds at the border of the event. Relatively small delays 
(Af_3 < 0.4) imply only a small reduction in the observed flux 
density (rff[Af_3] < 0.1), whereas Af_3 ^ 1 would be accompa- 
nied by significant or complete absorption of the signal. Inter- 
ferometric observations of the unpulsed continuum and obser- 
vations at other wavelengths will help to clarify the nature of 
the eclipses. 



a 
a 



26-Jan-2003 
UT 01:57 



0,1 0.2 0.3 0,4 

Orbital Phase 

Fig. 5. — Excess group delays of the signal of PSR J1701-3006B, measured 
on 2003 January 26 (UT time refers to orbital phase 0.25). The observation 
was centered at 1390 MHz with a bandwidth of 256 MHz and the data are pro- 
cessed in contiguous 360-s integrations. The error bars are twice the formal 
uncertainty in the pulse arrival times. The average value of the s/n within the 
eclipse region is 4.6 ± 0.6, whereas it is 5.7 ± 0.5 (1(7 uncertainty) outside. 

The occurrence of eclipses suggests that the orbital inclina- 
tion i is not small. For i = 60°, the median of all possible 
inclination angles, and an assumed pulsar mass of 1.40 Mq, 
Me,60 = 0.14 M0. For i > 30° the companion mass spans the 
interval 0.12-0.26 Mq, corresponding to a Roche lobe ra- 
dius in the range Rl = 0.26-0.34 Rq. Hence, independent of 

■^A possible exception is the pulsar C in the globular cluster M5 (Ransom, private 



the eclipse mechanism, the extension of the eclipsing cloud, 
^ 0.8 Rq, is larger than Ri and the cloud must be continu- 
ously refilled with matter released from the companion. The 
plasma density in the eclipse region is ~ 1.6 x 10"'^A?_3 
g cm"^ and, assuming isotropic emission, mass continuity im- 
plies that the donor star loses gas at a rate = AirR^pEVf ~ 
1 .0 X 10"'^ A?_3 V / 8 M0yr~^ , where v /,8 is the wind velocity at 
Re in units of 10^ cm s"^ (the order of magnitude of the escape 
velocity from the surface of the companion). 

If the companion is a helium white dwarf (whose maxi- 
mum radius is R^d = 0.04 R© for masses > 0.12 Mq and 
T < lO"' K, Driebe et al. 1998), and assuming isotropic emis- 
sion of the pulsar flux, a significant fraction / = (4% - 20%) x 
(3.7 X lO^'^erg s~^/E) of the energy deposited onto the compan- 
ion surface is necessary for releasing the observed Mc (where 
E is the spin-down power of the pulsar and 3.7 x 10^"* erg s~' 
its upper limit derived using the arguments of § 3). How- 
ever, the energy requirements are more easily satisfied for a 
non-degenerate bloated companion (as appears to be the case 
in most eclipsing binary pulsars, Applegate & Shaham 1994). 
For example, / = (0.04% -0.2%) x (3.7 x lO^-^erg s-^/E) for 
a donor with the radius of a main-sequence star of the same 
mass, that is, 3-10 times larger than that of a white dwarf. 
Mass loss from the donor star can be sustained by ablation of 
its loosely bound surface layers by the relativistic wind emit- 
ted by the pulsar This model has been successfully applied to 
explain the radio eclipses in close orbital systems having very 
light companions, e.g., the cases of PSRs B1957H-20 (Fruchter 
etal. 1990) and J205 1-0827 (Stappersetal. 2001). As with 
these other systems, the apparent mass-loss rate from the com- 
panion to PSR J1701-3006B is very small; the ablation time 
scale Tabi = X^cfio/Mc = xl40 Gyr, where x is the ionized frac- 
tion, is longer than the upper limit on the pulsar age (i.e. the 
cluster age) unless x < 0.09. 

Following an alternate interpretation, the PSR J1701-3006B 
system may more resemble that of PSR J1740-5340, where the 
effects of the pulsar irradiation are negUgible in triggering the 
eclipsing wind from the secondary star (D'Amico et al. 2001c) 
and the eclipses (or the excess propagation delays, sometimes 
seen far away the nominal phases of eclipse) are caused by mat- 
ter overflowing the Roche lobe of the donor star due to the nu- 
clear evolution of the companion (Ferraro et al. 2001). In that 
system, accretion of matter onto the neutron star is inhibited 
by the sweeping effect of the pulsar energetic wind, according 
to the so-called radio-ejection mechanism (Burderi, D'Antona 
& Burgay 2002). We note that J1701-3006B shares with PSR 
J1740-5340 (i) a companion significantly more massive than 
tiiose of PSRs B1957H-20 and J2051-0827, (ii) the occurrence 
of excess propagation delays at 1 .4 GHz which are much larger 
(up to ~ 1 ms vs few tens of /zs) than those observed in any 
of the systems having very low mass companions* and ( in) the 
presence of irregularities in the eclipses. 

A new class of eclipsing recycled pulsars having relatively 
massive companions (Me_60 = 0. 10-0.22 M0) is emerging from 
the globular cluster searches. Besides PSR J1701-3006B in 
NGC 6266 and PSR J1740-5340 in NGC 6397, there are PSR 
B1744-24A in Terzan 5 (Lyne et al. 1990), PSR J0024-7204W 
in 47 Tucanae (Camilo et. al. 2000) and PSR J2140-2310A in 
M30 (Ransom et. al. 2003a), whereas no similar system has 
been detected in the Galactic field to date. A simple explanation 
for the overabundance of evaporating "black widow" pulsars in 

communication) 
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globular clusters with respect to the galactic disk has been re- 
cently proposed by King, Davies & Beer (2003): namely, the 
current companion of most of the eclipsing pulsars in globulars 
would be the swelled descendent of a turn-off star which re- 
placed the original white dwarf companion of the pulsar in an 
exchange interaction in the cluster core. This scenario posits 
that the radio-ejection mechanism (Burderi et al. 2001) is now 
operating in all the eclipsing millisecond pulsars and provides 
an evolutionary basis for separating the systems with very low 
mass companion with respect to those having more massive 
donor stars; in the former, the mass loss would be driven by an- 
gular momentum loss through gravitational radiation, whereas 
in the latter the mass loss rate would be determined by the nu- 
clear evolution of the companion. 

The relatively massive systems in globular clusters are good 
candidates for optical detection of the donor star^ and follow-up 
observations. Unlike the Galactic eclipsing systems, their age, 
metallicity, extinction, distance and hence intrinsic luminosity 
and radius can be estimated from the parent cluster parameters 
(see, e.g., Edmonds et al. 2001a; Ferraro et al. 2001 ; Edmonds 
et al. 2002). In the case of J1740-5340 in NGC 6397, the com- 
panion is a red variable star of magnitude V ~ 16.5 (Ferraro 
et al. 2001) and stringent constraints have been set on the ef- 
fectiveness of the irradiation of the companion (Orosz & van 
Kerkwijk 2002), on the occurrence of the radio-ejection mech- 
anism (Sabbi et al. 2003) and on the evolutionary path of the 
system (e.g. Burderi, D' Antona & Burgay 2002; Grindlay et al. 
2002; Ergma & Sarna 2002). 

More recently, the companion of the millisecond pulsar 
J0024-7204W in 47 Tucanae has been optically identified with 
a blue variable star of mean magnitude V ~ 22.3, probably 
a heated main sequence star close to the center of the clus- 
ter (Edmonds et al. 2002). Unfortunately, the pulsar is weak 
and only occasionally detectable, which makes the system dif- 
ficult to characterize (Camilo et al. 2000). In the case of 
PSR B1744-24A in Terzan 5, the strong obscuration toward 
the Galactic center (> 7 mag in V) prevents detection of the op- 
tical counterpart, even with deep HST observations (Edmonds 
etal. 2001b). 

Consequently, PSR J1701-3006B is Ukely to be a primary 
candidate for improving the modeling of eclipsing millisec- 
ond pulsars with relatively massive companions. Indeed, PSR 
J1701-3006B seems to be a twin of PSR J0024-7204W in 47 
Tucanae, with similar orbital parameters and hence minimum 
companion mass (Camilo et al. 2000). Also the pulsar periods 
are comparable, 3.6 ms versus 2.4 ms. Moreover, unlike PSRs 
J1740-5340 and B1744-24A, both PSRs J1701-3006B and 
J0024-7204W reside well within one core radius of the parent 
cluster center and hence are in more similar environments. As- 
suming that the companion to PSR J1701-3006B has the same 
luminosity and colors as the companion to PSR J0024-7204W, 
its photometry would be feasible with deep exposures reaching 
V-magnitude 24.5. Photometry would of course be much easier 
if the companion fills its Roche lobe as is beUeved to be the case 
inPSRJ1740-5340. 

The X-ray counterparts of two of the five eclipsing MSPs 



with relatively massive companions (namely PSRs J1740-5340 
and J0024-7204W) have been identified using Chandra ob- 
servations (Grindlay et al. 2001a, 2001b). Their spectra ap- 
pear significantly harder than those of most other known X-ray 
counterparts to MSPs in globular clusters. That suggests (Ed- 
monds et al. 2002) that a non-thermal contribution to the X-ray 
emission, perhaps arising from shock interactions at the inter- 
face between the companion and pulsar winds, dominates over 
the thermal component seen in the other MSPs, which prob- 
ably originates from heated magnetic polar caps on the neu- 
tron star. The identification of the X-ray counterpart of PSR 
J1701-3006B and a comparison of its X-ray hardness ratio with 
that of the other MSPs in NGC 6266 would test the above pic- 
ture. Interestingly, a long Chandra pointing towards NGC 6266 
shows that it hosts the largest number of X-ray sources (with lu- 
minosity > 4 X 10^" erg s"^ in the 0.5-6.0 keV range) observed 
so far in a globular cluster (Pooley et al. 2003). Possibly among 
the 5 1 detected sources is a significant population of neutron 
stars, of which the 6 pulsars discovered so far are a manifesta- 
tion. 

7. CONCLUSION 

We have presented rotational and astrometric parameters of 
three binary millisecond pulsars located within 1 .8 core radii of 
the center of the globular cluster NGC 6266. One of these sys- 
tems, PSR J1701-3006B, displays ecUpses for ~ 20% of the 
orbit. In summary, we note that: 

1. The derived lower limits on the central mass-density 
(2.1 X 10^ MqPC"^) and the central mass-to-light ratio 
(M/C > 1.6 in solar units) of NGC 6266 are consistent 
with optical estimates. 

2. The large spread in the dispersion measures of the three 
millisecond pulsars is probably due to a significant gra- 
dient in the Galactic electron column density across dif- 
ferent lines-of-sight toward the cluster. 

3 . Even though the nature of the eclipses cannot yet be fuUy 

constrained, the relatively low mass-loss rate from the 
secondary star makes unlikely that PSR J1701-3006B 
will evaporate its companion. 

4. The lack of known isolated pulsars in NGC 6266 is un- 
Ukely to be due to chance or observational bias and sug- 
gests that the cluster is in a dynamical phase favoring 
formation over disruption of binary systems containing 
a neutron star. 
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'in fact, the optical identification of the secondary star has been recently reported also for two non eclipsing millisecond pulsars having companions with 
M^ eo ~ 0.2 Mq : PSR J0024-7204T in 47 Tucanae (Edmonds et al. 2003) and PSR 1911-5958A in NGC6752 (Ferraro et al. 2003; Bassa et al. 2003 ) 
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Table 1 

Observed and derived parameters for three pulsars in NGC 6266 



Pulsar name 


PSR J17()i-30()6A 


PSR J17()i-3006B 


P.SR.117()1-3()06C 


R. a. (J2000) 


17'" 01" 12?5127(3) 


01" 12?6704(4) 


l^h Qjm i2?8671(4) 


Decl. (J2000) 


-30° 06'30"13(3) 


-30° 06'49"04(4) 


-30° 06'59"44(4) 


Period, P (ms) 


5.2415662378289(16) 


3.5938522173305(14) 


3.8064243637728(18) 


Period derivative, P 


-1.3196(9) X 10*" 


-3.4978(7) X 10-1' 


-3.189(11) X 10-20 


Epoch (MJD) 


52050.0 


52050.0 


52050.0 


Dispersion measure, DM (cm"' pc) . . . 


115.03(4) 


113.44(4) 


114.56(7) 


Orbital period. Pi, (days) 


3.805948407(16) 


0.1445454304(3) 


0.2150000713(15) 


Projected semi-major axis, x (1-s) 


3.483724(8) 


0.252775(13) 


0.192880(12) 


Eccentricity? e 


<4x ur" 


< 7 X 10-5 


< 6 X 10-5 


Time of ascending node, Tasc (MJD) . . 


52048.5627980(15) 


52047.2581994(9) 


52049.855654(2) 


Span of timing data (MJD) 


51714-52773 


51714-52773 


51714-52773 


Number of TOAs 


80 


74 


73 


Rms timing residual (/is) 


21 


26 


32 


Flux density at 1400 MHz, S1400 (mjy) 


0.4(1) 


0.3(1) 


0.3(1) 


Derived parameters* 


Galactic longitude, / (deg) 


353.577 


353.573 


353.572 


Galactic latitude, h (deg) 


7.322 


7.319 


7.316 


Mass function, f], (Mq) 


0.00313392(2) 


0.00082999(13) 


0.00016667(3) 


Companion mass, Mc (Mq) 


> 0.20 


>0.12 


> 0.07 


Radio luminosity, L1400 (inJy kpc^) • • • 


19(7) 


14(6) 


14(6) 


Offset, ex (") 


19.2 


1.7 


10.5 



Note. — Figures in parentheses are twice the nominal tempo uncertainties in the least-significant digits quoted. 

^The 2(T upper limits on the orbital eccentricities were obtained using the tempo ELLl model, where Tasc and (eco%uj,e%iau>) are fitted (Lange et al. 
2001). The value given for PSR J1701-3006B is tentative as not all of the orbit is sampled. All the other parameters are derived using the standard (BT) 
binary model with e = 0. 

''The following formulae are used to derive parameters in the table: = {2tt / Pi,)^Tq = (Mc smi)^/(Mp +Mc)^, where T© = GMq/c^ = 4.925 /^s, Mp 
and Mc are the pulsar and companion masses, respectively, and / is the orbital inclination angle. Mc is obtained from the mass function, with Mp = 1.40 Mq 
(Thorsett & Chakrabarty 1999) and i < 90°. The assumed distance is that of the globular cluster, d = 6.9 kpc, and L1400 = Su(xid'^- 0± is the angular 
separation in the plane of the sky between the MSP and the center of NGC 6266 (Harris 1996, revision 2003). 



Table 2 

Encounter and disruption rates for binaries in four globular clusters 



Cluster 


Isolated PSRs 


Binary PSRs 




•?^dis.■ 


^form /^disr 


NGC 6266 





6 


1.4 


2.5 


0.57 


NGC 6752 


4 


1 


0.19 


3.4 


0.056 


M 15 


7 


1 


1.0 


5.5 


0.19 


47 Tucanae 


7 


15 


1.0 


1.0 


1.00 



Note. — T^-form is estimated as oc p^^r^, whereas 7?.disr oc Pq^'v' (see text for details). All the values are normal- 
ized to the parameters of 47 Tucanae. Central luminosity density po core radius are obtained from the catalog 
of Harris (1996, revision 2003). 



